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Abstract 

The most general QCD next-to-leading anomalous-dimension ma- 
trix of all four-fermion dimension-six AF = 2 operators is computed. 
The results of this calculation can be used in many phenomenologi- 
cal applications, among which the most important are those related 
to theoretical predictions of K°-K° and B°-B° mixing in several ex- 
tensions of the Standard Model (supersymmetry, left-right symmetric 
models, multi-Higgs models, etc.), to estimates the B®-B® width dif- 
ference, and to the calculation of the 0(l/mf ) corrections for inclusive 
6-hadron decay rates. 



1 Introduction 



Theoretical predictions of several measurable quantities, which are relevant 
in K-, D- and 5-meson phenomenology, depend crucially on the matrix 
elements of some AF = 2 four-fermion operators. Examples are given by 
FCNC effects in SUSY extensions of the Standard Model (or other 

models such as left-right symmetric || or multi-Higgs models), by the B®- 
width difference [|J, and by the 0(1 /ml) corrections in inclusive b-hadron 
decay rates (which actually depend on the matrix elements of several four- 
fermion AF = operators ||). In all these cases, the relevant operators 
have the form 

Q = C a ^(b a Tq /3 )(b p Tq a ), (1) 

where T is a generic Dirac matrix acting on (implicit) spinor indices; a-a are 
colour indices and C al3pa is either 5 al3 5 pa or 5 aa S pl3 (for the 1/m 3 corrections 
to the inclusive decay rates the flavour structure has the form (bq)(qb))!^ 

All the operators discussed in this paper appear in some "effective" the- 
ory, obtained by using the Operator Product Expansion (OPE). As a conse- 
quence, in all cases, three steps are necessary for obtaining physical ampli- 
tudes from their matrix elements: 

i) matching of the original theory to the effective one at some large energy 
scale; 

ii) renormalization-group evolution from the large energy scale to a low 
scale suitable for the calculation of the hadronic matrix elements (typically 
1-5 GeV); 

iii) non-perturbative calculation of the hadronic matrix elements. 

In this paper we present a calculation of the two-loop anomalous dimen- 
sion matrix relevant for AF = 2 transition amplitudes. This matrix can be 
used for the Next-to-Leading Order (NLO) renormalization-group evolution 
of the Wilson coefficient functions of the effective theory from the large to 
the small energy scale, step ii). The anomalous dimension matrix includes 
leading and sub- leading corrections of order a s and a 2 s . The calculation has 
been performed in naive dimensional regularization (NDR) and we give re- 

1 All the formulae of this paper refer to the AB — 2 case. Their extension to generic 
AF = 2 transitions is straightforward. 
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suits for different renormalization schemes. We have also verified that the 
results obtained in different schemes are compatible. We give many details 
on the calculation itself, on the definition of the renormalized operators, on 
the relation between different renormalization schemes (and on the role of 
the corresponding counterterms) , on the gauge invariance of the final results 
etc. We also present a list containing the contribution of all the Feynman di- 
agrams to the anomalous dimension matrix. The list may be useful to check 
our results and for further applications. In this paper we have preferred to 
give the results for the one- and two-loop anomalous- dimension matrix with 
as many details as possible and postpone the phenomenological applications 
of the results given here to further publications. In section |5|, the reader who 
is not interested in the theoretical and technical aspects of the calculations 
can find the final results for the anomalous dimension matrix of the operator 
basis defined in eq. ( |T3| ) of subsec. . 

Besides presenting the results for the AB = 2 (AF = 2) operators dis- 
cussed here, we take the opportunity to clarify several issues related to the 
regularization and renormalization dependence of the operators and of the 
corresponding Wilson coefficients. In particular, we discuss in detail the 
problems related to the precise definition of the so-called "MS schemes", 
which, for composite operators, are not uniquely defined, even for a given 
regularization 0-[Tj|. We also examine the equivalence, and differences, 
of the most popular renormalization schemes, and the subtleties related to 
Regularization- Independent renormalization schemes (RI) |14{| . 

The paper is organized as follows. In sec. [| we introduce the operators 
relevant for physical applications and define the operator basis for which the 
anomalous dimension matrix will be given; a general discussion on the Wilson 
coefficients, their scheme- dependence and renormalization-group evolution 
will be presented in sec. |3|; the strategy for the calculation of the anomalous 
dimension matrix in the MS and RI schemes, to be defined below, is given 
in sec. the final results, together with the one-loop matrices necessary 
to change renormalization scheme, are also given for some relevant cases in 
sec. H. 
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2 Four-fermion operators 



We start this section by introducing the operators that we have in mind in 
view of future applications; we then illustrate the chiral and Fierz properties 
of the relevant operators which are used to derive the general form of the 
mixing matrix. 



2.1 Operators relevant for physical applications 

In this subsection, we present a list of operators which enter the calculation 
of the physical quantities mentioned in the introduction |: 

1) FCNC in SUSY extensions of the Standard Model: 

For K- and .B-meson transitions, these effects have been recently an- 
alyzed in detail in a series of papers, see for example refs. The 
relevant operators which enter the effective Hamiltonian are 

Q x = b a lfl (l- l5 )q«b\(l- l5 )qP, 

Q 2 = - b a (l- l5 )q a bP{l- l5 )qP, 

Q 3 = b a (l- l5 )q^¥{l- l5 )q a , (2) 

Q 4 = ^(1-^)^(1 + 75)/, 

Q 5 = 6 Q (i-7 5 )g /3 & /3 (i + 7 5 )g a , 

together with the operators <3i,2,3 which can be obtained from the op- 
erators <3i,2,3 by the exchange (1 — 75) (1 + 75). 

2) The B s -B s width difference Ar# s : 

At lowest order in l/m&, by using the OPE, the width difference Ar# s 
can be written in terms of two AB = 2 operators 0] 

Q = 67^(1-75)567^(1-75)3, 
Q s = 6(l- 7 s)s6(l- 7B ) S . (3) 

where, since the fermion bilinears are colour singlets (67^(1 — 75 )s = 
6 a 7 A1 (l — 75)5"), the colour indices have not been shown explicitly. 



Here and in the following we adopt the same notation as in the original papers. 
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3) Heavy-hadrons lifetimes (tb, tb s , taJ: 



In this case, the 1 / ml corrections to the lifetime, due to Pauli interfer- 
ence and W -exchange, can be written in terms of four operators || 

6T7*(l-75)g 97 M (1 -75)6, 
6(1-75)99(1 + 75)6, (4) 
U A lfl (l- l5 )qqt A Y(l-l 5 )b, 
U A (l- l5 )qqt A (l + l5 )b. 

where an implicit sum over colour indices is understood. The operators 
above are AB = operators. They contribute to the decay rates of the 
.B-mesons (and A^s) not only through the so-called "eight" diagrams, 
but also through tadpole diagrams, in which the light- or heavy-quark 
fields are contracted in a loop. These "non-spectator" diagrams mix 
the operators of the basis @ with the lower dimension operators 66, 
bD 2 b, ba^G^b, qq, etc. The mixing matrix is, however, triangular. 
Thus, it is possible to compute separately, at the NLO, the 4x4 sub- 
matrix related to the mixing of the operators appearing in (^) among 
themselves. For this sub-matrix, the Feynman diagrams entering the 
calculation are the same as those relevant for the AB = 2 operators. 



n q 

U V-A 
U S-P 

rpq 

2 V-A 

rpq 

1 S-P 



2.2 Chiral and Fierz properties of the operators 

The operators considered in l)-2) can be expressed in terms of linear com- 
binations of independent operators, defined by their colour-Dirac structure, 
belonging to some basis. In case 3), for the sub-matrix considered here, 
the same colour-Dirac structure (with obvious replacement of the flavour in- 
dices) can also be used. The choice of the basis of reference is, however, 
arbitrary, and different equivalent possibilities exist. We first present one 
possible choice, which we find particularly convenient to discuss the chiral 
properties of the operators: 

Qv L v L = ^(1-75)^2^(1-75)^ 

QvlVr = ^(1-75)^2^.(1+75)^ 

Qrl = ^(1+75)^2^(1-75)^ 
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Qll = ^(1-75)^2^(1-75)^ 

QT L T L = ^i<v(l ~ 75)^2 ^3^(1 - 75)^. 



Qv.v, = ^?7m(1 - 75)^ *7m(1" 75)^4 (5) 

QvlVr = ^(1-75)^1^(1 + 75)^4 

Q«l = ^(1 + 75)^^(1-75)^ 

Qll = ^(1-75)^^(1-75)^ 

Qt £ T £ = V>?0>(1 - 75)^2 ^3^(1 - 75)^4 , 

where cr^ = 1/2 [7^,7^]. In (0), the flavours ipi-ipi are all different and the 
operators belong to irreducible representations of the chiral group. To these 
10 operators, we have to add those which can be obtained by exchanging left- 
with right-handed fields. Since, however, strong interactions cannot change 
chirality Q, the second set of operators does not mix with the operators defined 
in (H) and, because parity is conserved, the Anomalous Dimension Matrix 
(ADM) is the same in the two cases. Thus, in the following, we will only 
consider the operators of eq. @. 

The previous considerations hold only if one uses a renormalization pre- 
scription which preserves chirality (in this respect parity is never a serious 
issue). It often happens, e.g. in dimensional schemes such as the t'Hooft- 
Veltman MS one (HV), that the renormalization procedure violates either 
chirality, or (and) other symmetries that are manifest at the tree level, for 
example the Fierz transformation properties. In order to simplify the pre- 
sentation of the results, we will use in the following a renormalization scheme 
which preserves all the relevant symmetries (chirality and Fierz). With such 
a choice, it is then sufficient to consider the basis (|5|). This renormaliza- 
tion scheme has been recently called the Regularization Independent (RI) 
scheme |H| (MOM in the early literature) to emphasize that the renormal- 
ization conditions are independent of the regularization, although they de- 
pend on the external states used in the renormalization procedure and on 
the gauge. The RI scheme offers also a great computational advantage in 
the calculation of the counterterms which contribute at the two-loop level: 



as demonstrated in subsec. |4.2|, in this scheme it is not necessary to iden- 

3 In mass-independent renormalization schemes, such as the RI schemes discussed in 
this paper, chiral-symmetry relations, which can be derived in the massless theory, remain 
true also in the massive case. 
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tify and subtract separately the counterterms relative to the mixing with 
the "Effervescent Operators" (EOs) which appear in dimensional regular- 
ization 0— [13]. The relation between the operators renormalized in some 



RI scheme and, for instance, those of the standard MS schemes can then 
be easily found with a simple one-loop calculation. Finally, the RI scheme 
allows to use, without any further perturbative calculation, the matrix el- 
ements of the operators computed in lattice simulations and renormalized 
non-perturbatively |15|-p!7|. 



Chiral symmetry, and Fierz rearrangement, have further consequences, 
since they forbid the mixing between some of the operators appearing in 
eq. For this reason the ADM, 7, is a block-matrix which only allows 
mixing between sub-sets of the possible operators. In the convenient rep- 
resentation in which the operators appearing in (|5|) are components of row 
vectors 

I) Qi = (Qv L v L , Qv l v l ), 
II) Qn = (Qv L v R ,Qv L v R ), 

III) Qui = (Qrl, Qrl), 

IV) Qiv = (Qll,Qll,Qt l t l ,Qt l t l ), 

Fierz rearrangement imposes the following restrictions on the form of the 
mixing matrix: 

• For set I) the structure is given by 

• If the mixing-matrix for II) is given by 

- — ( An B n \ , , 

1U =\c n D n )> ^ 

then we have 

7m = o a.). (8) 
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In order to discuss Fierz rearrangement for the sector IV), we introduce 
the Fierz transformation matrix for the sub-basis Qiv 



( -\ 

\ o I 





V 6 



6 




(9) 



I 0/ 

The anomalous dimension matrix must satisfy the relation 

j iv = Flw? 

The mixing matrix has a simple form in the Dirac-Fierz basis 



(10) 



Qf 


— (Qi i Q2 1 Qi -i Q2 


Qf 


= Qll + -^Qt l t l , 


Ql 


= Qll - -^Qt l t l , 


Qf 


= Qll + -^Qt l t l , 


Ql 


= Qll - t^Qt l t l ■ 



(11) 



In this basis, the ADM can be written as 



liv 



( 


Aiv 


Biv 


Civ 


Div 


\ 




Eiy 


Fry 


Giv 


Hiv 




Civ 


-Div 


Aiv 


—Biv 




V 


—Giv 


a iv 


—Eiv 


Fiv 


J 



(12) 



In summary, we have seen that the ADM for all the operators appearing 
in (|) can be expressed in terms of 14 quantities, i.e. Ai, Bi, An, . . ., Hiv- 



2.3 The Fierz basis 

In the case in which ipi = ip 3 (or ip 2 — ^4)5 n °t a U the operators in (||) are 
independent. In order to take into account the simplifications occurring in 
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this particular case, it is more convenient to give the results in the (Fierz) 
basis 



Qt 
Qt 
Qt 



Qt 



(Qt-, Qti Qti Qai Qt) 



Q 



X - (^(1 - 75)^2 $7m(1 + 75)^ ± (V>2 <- ^4 

i + 75)^2 *(1 - 7s)< ± (^2 <- V4)) 

I (^(i _ 75 )^ ^(1 _ 75 )^ ± <_> 

^ (^^VC 1 - 75)^2 V>3<V( 



(13) 



75)^4 ± (^2 <"> ^4 



In this case, the operators do not belong, in general, to irreducible repre- 
sentations of the chiral group, e.g. both right- and left-handed ip2 fields 
appear in Qf . The AB = 2 operators are obtained from the Qfs, by taking 
ipi = V's = b and ip2 = ip4 = Q (with this choice of flavours the Q^ vanish). 
In the basis ([T3]), the ADM has the form 



7± = 



( A± 














\ 





B 


±c 













±D 


E 



















F ± 


G ± 




V 








H± 


J± 


/ 



(14) 



and there is no mixing between the Qf and the Q { operators. 



The correspondence between the operators of the basis (0) and the op- 
erators which are relevant for the physical applications listed in subsec. ^[T] 
is the following 



Qi Qt > Q2 
Qa — > Qt 5 Q5 — > 

Q Qt > Qs - 



* Qt 1 

I — 
"2 

Qi > 



2 1 



(15) 
(16) 
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3) Oy_ A — > Qi + Q 1 , Q q s-p — » Qt + Q3 > 

n-^K'-wjrt-kk + w)* (17) 

TL P - (Qi + g 3 " ) - \ [Qi - Q2) ■ 
3 The Wilson coefficients 

The general method for the calculation of the Wilson coefficients, and a de- 
tailed discussions on their renormalization-scheme dependence, can be found 
in the literature [|l4[]-[|n|. In this section, we only summarize the main formu- 
lae which are necessary to present our results. We also take the opportunity 
to clarify some important subtleties about the renormalization-scheme de- 
pendence. 



3.1 Effective theories and Wilson coefficients 

In all cases of interest, the matrix elements of the effective Hamiltonian can 
be written as 

(F\H eff \I) = Y,(F\Qi(v)\I)C^), (18) 

where the Qi(fi)s are the relevant operators renormalized at the scale \x and 
the Cj(//)s are the corresponding Wilson coefficients. We represent the op- 
erators as row vectors Q, as in subsec. ^2| , and the coefficients, C(fi), as 
column ones. The vectors C(/x) are expressed in terms of their counter-part, 
computed at a large scale M, through the renormalization-group evolution 
matrix W\pt, M] 

C(fi) = W[fx,M}C(M). (19) 

The initial conditions for the evolution equations, C(M), are obtained by 
matching the full theory, which includes propagating heavy-vector bosons 
(W and Z°), the top quark, SUSY particles, etc., to the effective theory 
where the W, Z°, the top quark and all the heavy particles have been re- 
moved simultaneously. In general, C(M) depend on the definition of the 
operators in a given renormalization scheme. The coefficients C(fi) obey the 
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renormalization-group equations: 



9 Of 



_d_ 

da. 



+ /3 A (« S )A--^— 



C(t,a.(t),A(t)) = 



(20) 



where t = ln(M 2 //i 2 ). The term proportional to (3\, the /3-function of the 
gauge parameter \(t) (for covariant gauges), takes into account the gauge 
dependence of the Wilson coefficients in gauge- dependent renormalization 
schemes, such as the RI scheme [14[ || []. This term, the role of which 
will be discussed extensively in sec. |], is absent in standard MS schemes, 
independently of the regularization which is adopted (NDR, HV or DRED 
for example) [f7[-||13||. The factor of 2 in eq. (|20| ) normalizes the anomalous 



dimension matrix as in refs. [fl4|-[113|. To simplify the discussion, we only 
consider the case where there is no crossing of a quark threshold when going 
from M to fi. The relevant formulae for the general case can be found in 



refs. Q-|3|. 

At the next-to-leading order, we can write 

W[pi,M] = M[^]U[fi, MjM'^M] 

where U is the leading-order evolution matrix 

7 (0)T /2/3o 



and the NLO matrix is given by 



a s (M) 



(21) 



(22) 



M[fx] = 1 + J[A(a*)] 



(23) 



By substituting the expression of the C (/x) given in eq. ( |19|) in the renor- 
malization-group equations (|20|), and using W[/i, M] written as in eqs. (pT])- 



3[), we find that the matrix J satisfies the equation 



j + ^x dj 



J 



7 



(0)T 



2/?o 



A_.(0)T _ 7 



(1)T 



2Po 



(24) 



In the following, we will denote by A = 1 the Feynman gauge and A = the Landau 



gauge. 
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In eqs. (p^) and fl2"4|), /5o, /3i and are the first coefficients of the (3- 
functions of a s and of A, respectively; 7^ and 7^ are the LO and NLO 
anomalous dimension matrices to be defined in sec. U is determined by 
the LO anomalous dimension matrix 7^°) and is therefore regularization and 
renormalization-scheme independent; at this order, XdJ/dX is also regular- 
ization (but not renormalization) scheme independent; the two- loop anoma- 
lous dimension matrix 7W, and consequently J and W\p,, M], are, instead, 
renormalization-scheme dependent . 



3.2 Coefficient functions and scheme dependence 

In this subsection, we recall some basic aspects of the calculation of the 
Wilson coefficients and discuss in detail the issues of the regularization and 
renormalization dependence of the coefficients and of the corresponding op- 
erators. We believe that this discussion may be useful to clarify some mis- 
understandings that can be found in the literature. 

In order to compute the Wilson coefficients at a large energy scale fi ~ M, 
we should consider the full set of current-current, box and penguin diagrams 
in the full theory, i.e. with propagating heavy particles, including the 0(a s ) 
corrections. To date, for the AF = 2 transitions, this part of the calculation 
has been carried out only in the Standard Model and 2HDM cases H-|ITJ . 



In the full theory, the direct calculation of the current-current, box and 
penguin diagrams, including 0(a s ) corrections has the form 

<«.//> ~ (Q i0)T ) ■ [T {0) + = <<f (/,)> ■ , (25) 

where (Q^ T ) are the tree-level matrix elements and the vector T^> depends 
on the external quark (and gluon) states chosen for the calculation. By 
inserting the renormalized operators of the effective Hamiltonian, we then 
compute, at order a s , the one-loop diagrams between the same external states 
as in the full theory, using the same regularization. In this case we obtain f] 

a, „\ , - (0)v 



(Q(/i)) = (l + ^r)<QW>. (26) 



5 The most convenient method to define the matrix elements is by projectors on the tree- 
level colour-Dirac structures of the operators belonging to the four-dimensional basis |16| . 
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C(M) = f (°) + ^ (f « - f T f (°>) . (27) 



The coefficients C(/x) are obtained by comparing eq. (p5|) with eq. (|26|); if 
we (formally) choose the renormalization scale fj, = M, all the logarithms 
related to anomalous dimensions of the operators disappear and 

and r T depend on the external states. However their difference depends 
only on the renormalization scheme, but not on the external states. For this 
reason, the dependence on the external momenta (~ ln(— p 2 )) of f in C(M) 
and cancels out, for details see refs. [11]— [13]. In the following f and 
denote only non-logarithmic terms. 

For given external states, and for a given gauge, the matrix f completely 
specifies the renormalization scheme (MS, RI, etc.). In this respect, all the 
renormalization schemes, including the MS ones, are regularization indepen- 
dent. The MS schemes simply amount to some specific choice of f. This 
is also demonstrated by the following observation: even when the regular- 
ization is specified, for example the NDR one, the so-called "MS scheme" 
is not unique. The renormalized operators, and consequently f, depend in 
general on the basis chosen in the regularized theory to implement the min- 
imal subtraction procedure, the projectors, i.e. the definition of the EOs, 
etc. [|TT[] — pT3f] . Thus, in order to define completely the "MS scheme", we 



should specify all the variables (regularized basis, EOs, etc.) entering the 
calculation. In practice, this is equivalent to fix f, i.e. the renormalization 
prescription. Summarizing, the regularization dependence must always be 
understood as a "renormalization-scheme dependence" . In all the renormal- 
ization schemes (both MS and RI), the same information is contained in f, 
once that the external states and the gauge are specified: we will then use 
eq. fl26|) to define the renormalized operators. The explicit expressions of the 
matrix f in the different schemes (and the corresponding states and gauge) 
can be found in sec. 01. 



In subsection J4.2L it will be shown that the combination 



G = 7« - \r, 7 (0) 1 - 2/3f - 2(3° X X^- (28) 

- J OA 

is renormalization-scheme independent. It can be easily shown that a conse- 
quence of eq. (|28f) is the independence of the combination 

Jri = J + r T (29) 
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of the renormalization scheme (but not of the external states and, in general, 
of the gauge on which f is computed), see also refs. [14]— [13] . 

The independence of Jrj (and of G) of the renormalization scheme means 
the following. As mentioned above, for given external states, and in a given 
gauge, we compute the matrix element of the renormalized operators. These 
operators are renormalized in some scheme, for example one of the possible 
MS schemes. If we change scheme, f and J will accordingly change, whilst 
Jri will remain the same. Using eqs. ( p4| ) and (|28|), we find indeed 



J, 



RI 



A) d\ 



J, 



7 



(0)T 



RI, 



2 A) 



2pf 



(0)T 



Wo 



(30) 



The scheme independence of G and 7^°- ) guarantees the independence of the 
solution of eq. (p0[). In turn, this implies the renormalization-scheme inde- 
pendence of the matrix J RI . Note also that, since Jj^g is gauge independent, 

dj RI /d\ = dfyd\. 

The renormalization-invariant properties discussed above can be used to 
introduce schemes which respect all the symmetries of the tree-level theory. 
Using eqs. (pT|) (p3|) and (|27f) , we introduce a new set of Wilson coefficients 
C\M) 



C(n) = M[/4 U[n,M] N-^M] C'(M) 
where M\p] has been defined in eq. (^3|) and 

a s (M) 



N[M] = 1 + 



47T 



J + f 1 



C'{M) 



f (0) + ^±f(i) 
47T 



In the above equations we have neglected higher order terms in a s 
suitable change of the renormalization scheme, corresponding to 



(31) 

(32) 
By a 

(33) 
(34) 



Equation (|34f) has the following interpretation: it corresponds to the general 
expression (pl|), with the matrices M[fi] and M[M] given in terms of Jri, 



V T (») = Q r (/i) 1 - 



Att 



and — > N[/j] one gets 

C\n) = N[/j]U[h, MjN-^MjC'iM) . 
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which satisfies eq. (|30|). In this renormalization scheme, G = 7^, since f = 
Fri = 0. Clearly, the result is independent of the renormalization scheme of 
the original operators Q(fi)- The scheme dependence is implicitly contained 
in the states and the gauge on which f is computed. Thus, in the following, 
we will call FRI (LRI) the scheme with matrix elements computed in the 
Feynman (Landau) gauge. 

A remark is in order at this point. In refs. [I8|— [20|, for AB = 1 transi- 



tions, the authors used the so-called regularization-scheme-independent co- 
efficients (corresponding to our C"(/x)) introduced in ref. [fy| | and computed 
in ref. [[Tj|. In this reference, the preference for using this particular renor- 
malization scheme was justified with the argument that also the operator 



matrix elements computed with factorization (and used in [18]) are scheme 
independent. It was therefore argued that the regularization- independent 
coefficients are more suited to obtain the physical amplitudes. This argu- 
ment is clearly illusory: the coefficients, though regularization independent, 
depend on the external states and on the gauge at which the renormaliza- 
tion conditions have been imposed. There is no way to match the external 
quark and gluon states, used in the perturbative calculation of the Wilson 
coefficients, to the hadronic states on which the operator matrix elements 
are computed (not to speak about gauge invariance). A similar argument 
applies to the scheme used in ref. ||21[| , where the authors try to get rid of 
the \x dependence of the non-leptonic amplitudes computed with factoriza- 
tion. In any case, in the absence of a consistent calculation in which both 
the coefficients and the matrix elements of the operators are computed with 
the same renormalization, a preferred scheme does not exist. 



4 Anomalous dimensions at one and two loops 

In this section we recall the procedure for the calculation of the anomalous 
dimension matrix 7 in dimensional regularization. The section is divided in 
two parts: in the first part, we introduce the general formulae which define 
the anomalous dimension matrix, in the second we give the practical recipe 

6 Indeed the renormalization scheme of ref. [[ll] has never been completely specified, 
because the external states, on which the renormalization conditions were imposed, have 
not been given explicitly. 
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to compute it at the NLO. 



4.1 General definitions and scheme dependence 



The ADM for the operators appearing in the effective theory is given by the 
matrix 

7 = 2Z~V 2 :^, (35) 



where Z = Z(a s , A), defined by the relation 



Q 



Q 



(36) 



gives the renormalized operators in terms of the bare ones. A is the renormal- 
ized gauge parameter on which, in general, Z may depend, see for example 
ref. 0. Note that in refs. the dependence on the gauge parameter was 

ignored because, in MS schemes, Z is gauge independent. In this work, since 
we will compare anomalous-dimension matrices between schemes in which Z 
can be gauge dependent, such as the RI scheme, this dependence has to be 
taken explicitly into account. 



In dimensional regularization, using eq. 

7 = 2Z- 1 



one gets 



d . d . 

(-ea s + (3{a s )) - — Z + \(3 x {a s ) — Z 
oa. OA 



(37) 



where e = (4 — D)/2. (3{a s ) and (3\{a s ) are the (3 functions which gov- 
ern the evolution of the effective coupling constant and renormalized gauge 
parameter A, respectively 



>da s 
djj? 



with 

(5{a s 



'Air 



a" 



(4tt) : 



(5(a s 



+ 0{a 



4 ) 

si 1 



dX 
d\i 2 



%0l + O{a 



(38) 



(39) 



Po, Pi and fix are given by 



Pi 



(UN -2nj 
3 

N /13 



01 



34 Ar2 10 Ar iV 2 -l 

— Nn f - -n f , 

3 3 N 



(40) 
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where rif is the number of active flavours. The strong coupling constant 
a s and the gauge parameter A are renormalized in the MS scheme. This 
does not imply that also the operators must be computed in MS, because 
the definition of the composite operators is an independent step that has 
nothing to do with the procedure which renormalizes the parameters of the 
strong-interaction Lagrangian. 

From eq. ([37]), by writing 7 and Z as series in the strong coupling constant 
7 = ?7 (0) + 7^%7 (1) + ---, (41) 



An (An] 



and 



^=1 + ^ (1) + 7^U (2) + ---, (42) 
An (An) 2 

we derive the following relations 

7<°> = -2eZ« (43) 

and 

T (D = _ 4e ^(2) _ 2 /3 Z^ + 2eZ^Z^ - 2(3 Z^ - 2/3 A °A^— . (44) 

OA 

We can expand in eqs. ( f43~|) and (f44[ ) in inverse powers of e 

Z« = £ (-Y • (45) 

The requirement that anomalous dimension is finite as e — * implies a 
relation between the one- and two-loop coefficients of Z (note that in all the 
regularizations z[^ is gauge invariant for gauge invariant operators) 

AZi 2) + 2[3 z[ 1) - 2Z[ 1) z[ 1) = , (46) 

which can be used as a check of the calculations. In addition, from the 
eqs. (filf ) and (|44]) we obtain 



7 (°) = -2Z{ 1] (47) 
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and 



7 (1) = "4^ 2) - 2/3 ^ 1) + 2(Z< 1) Z< 1 > + Z^Z?) - 2#A^ • (48) 

Thus, it is sufficient to compute the pole and finite part of Z^ and the 
single pole of Z@\ together with (3q and /3°, in order to obtain the two-loop 
anomalous dimension. Note that the last term in eq. ([48]) is absent in refs. J?j- 
[ |T3|j . Eq. ( f4*8|) tells us how to derive In dimensional regularizations, such 
as HV, NDR or DRED, the calculation is complicated by the presence of the 
so-called "effervescent" operators, which appear in the intermediate steps of 
the calculation [|, 0]. The EOs are independent operators which are present 
in D dimensions but disappear in the physical basis of the 4-dimensional 
operators. Because of the presence of the EOs, the products of the matrices 
Zj in eq. fl48|) have to be done by summing indices over the full set of 
operators, including the EOs. Only at the end of the calculation we can 
restrict the set of operators to thoses of the physical 4-dimensional basis. 
As explained below, the identification of the EOs, and of the corresponding 
mixing matrix, can be completely avoided in RI schemes. 



4.2 Extraction of the the one- and two-loop anomalous 
dimension matrix 

We now derive the general expression of the coefficients Zj s in an arbitrary 
renormalization scheme, as obtained by using dimensional regularization. 
The derivation is general and, with trivial modifications, holds also with 
other regularizations, such as for example the lattice one |14] . Let us consider 



the matrix elements of generic bare operators, denoted as Qb, computed in a 
covariant gauge, between assigned quark and gluon external states. We define 
the matrix elements of Qb as the 1PI bare Green functions Tq b multiplied 
by the renormalization constants of the external fields 

{Q B ) = Z^T Qb , (49) 

where Z^ will be defined below. By calling «o the dimensionless bare coupling 
constant, for p 2 = — fi 2 , where p 2 denotes generically the squared momentum 
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of the external states, we have [] 



(Qb) 



(Q {0) ) ■ (50) 



(Q^) are the tree- level matrix elements of all the operators of the regularized 
theory, including the EOs. At one loop, for gauge-invariant operators, the 
matrix A\ is gauge and regularization independent, whilst Aq can be written 
in the form 

8 A 

io(Ao) = io(0) + A ^ (51) 

where A is the bare gauge-parameter. Note that also dA /d\ is regulariza- 
tion independent. 

In eq. (|5UD, we substitute the bare parameters a and A with their 
renormalized counter-parts, according to 

a. = a(i-£? + ...), - = -(i-St + -)" < 52 > 

At the NLO, and taking into account that A% is gauge invariant, we can 
ignore all other terms which relate the bare and the renormalized coupling 
constant and gauge parameter. 

For a given, generic renormalization scheme, we can write the following 
relation between matrix elements 

(Q R ) = Z-\Q B )= (1 + ^(3(0)), (53) 

where the matrix f defines the renormalization scheme. With a little algebra, 
a comparison of eqs. (|50D, fl5"2"D and (|53|) gives the mixing matrix Z in terms 

7 When working in the MS scheme, it is convenient to express the poles in terms of 
1/e = 1/e — 7_e + ln(47r), where 7^; is the Euler gamma. The formulae below are valid also 
in the MS scheme if one interprets 1/e as 1/e. 
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of An, . . . B2 (we list only the terms which are needed at the NLO) 



Z[ 2) 



An -r 



Z? = A, 



B, - A 1 f - P A - (3° X X 
B 2 -f3 A 1 . 



dA, 
OX 



(54) 



7W is then readily obtained by substituting the relations ( |5"4] ) in eq. (£8). To 
this purpose, we express the regularization and renormalization-independent 
combination (obviously y^ = —2Ai) 



G = 7 (1) 

in terms of the matrices A , 

1 



r,7 
B 2 



(0) 



(JT 

2/3 r - 2Pl\— 



G 



B 1 -- (iiio + AoA,) - ^[3 A - \p°x^ 



(55) 



(56) 



Equation (|56D demonstrates that G is renormalization-scheme independent 
since the r.h.s. does not depend on f. G is also regularization independent 
as the following, simple argument demonstrates. Let us compute the matrix 
elements of the renormalized operators using two different regularizations, 
but the same external quark and gluon states and in the same gauge. Irre- 
spectively of the regularization used in the calculations, we can define the 
operators in the same renormalization scheme in terms of the mixing matrix 
f in eq. (|53|). Since the renormalized operators are the same, they obey the 
same renormalization-group equations. Thus, not only f, but also is the 
same in the two cases. This demonstrates that the l.h.s. of eq. fl56|) is also 
regularization independent. 

Note that the last term of eq. Q55D is absent in refs. This is 

because in all MS schemes the derivative df/dX is regularization invariant, i.e. 
it is the same for two different MS regularizations. Thus, in these schemes, 
the difference between the two-loop ADMs is given by: 



A7 



(i) 



Af,7 (0) +2/3 Af 



(57) 



On the other hand, the combination G depends on the external states used 
in the calculation, and on the gauge, because f depends on these variables. 
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The RI scheme is defined, for given external states and at a fixed gauge, by 
the condition f = 0. Thus, in this scheme, G coincides with the two-loop 
anomalous dimension. 

Equation (|56|) provides also a practical method to compute the ADM in 
terms of the one- loop matrices A\ and Aq and of the two-loop pole term B\. 
In refs. and [fL3| , it was demonstrated that the equation which allows to 



compute in terms of the one- and two-loop renormalization matrices is 



valid diagram-by-diagram. In ref. [O] it was also shown that the relations be- 



tween the anomalous dimensions in different regularizations/renormalizations 
can also be established on a diagram-by-diagram basis. Using these obser- 
vations, and eq. (|56|), we give the recipe to obtain in a very simple way the 
anomalous dimension matrix in the RI scheme, 7^: 

i) Choose the set of external states, and the gauge, which define the RI 
scheme of interest. 

ii) Compute a given two-loop diagram where the bare operator is inserted. 

iii) Subtract to it the result obtained by substituting, to any internal subdi- 
agram, one half of the amplitude of the corresponding one-loop diagram 
computed at p 2 = — /i 2 : 

-W 0) >> (58) 




i.e. one half of the contribution of the subdiagram to A and A±. 

iv) When the internal subdiagram contributes to the renormalization of a s 
and of the gauge parameter A, apply iii) by inserting, in the two-loop 
diagram, only the divergent part of one-loop diagram (always with a 
factor 1/2). This rule corresponds to the choice of renormalized a s and 
A in the MS scheme. 

v) The coefficient of the single pole obtained from steps i) - iv) is the 
contribution of the given two-loop diagram to the combination (|56|) . 



With this procedure, we do not need to isolate the EOs from the operators 
of the 4-dimensional basis. The reason is two-fold. On the one hand, in the 
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RI scheme, the counterterms corresponding to the EOs and to the operators 
of the 4-dimensional basis (i.e. the combination AiAq + AqA\ of eq. (|56[) ) 
are both subtracted with the same factor 1/2. As shown in eq. (plj ) below, 
in the MS scheme, instead, different factors, namely 1 and 1/2, enter the 
subtraction for the 4-dimensional and effervescent counterterms. Moreover, 
the subtracted diagrams, as obtained from steps i) to v), only contain simple 
poles or finite terms, while the double poles completely cancel out. Thus, 
the projection on the physical 4-dimensional basis cannot give rise to further 
single-pole terms due to the EOs. 

For completeness, we now give the definition of the quark wave-function 
renormalization which we used in the RI scheme. We introduce the two-point 
Green function, computed in the same gauge as the RI scheme at hand, 

where the trace is taken over colour and spin indices. The renormalization 
condition for the quark fields is given by 

V)!V(p 2 = /i 2 ) = i ■ (60) 

Note that, in the calculation of the four-point Green functions, in any given 
scheme, different choices of the wave-function renormalization correspond 
to different choices of the quark external states. Thus, they also imply, in 
practice, different definitions of the renormalized operators. Obviously, all 
these choices are equivalent in principle, and they do not affect the calculation 
of the physical quantities at the order we are working. In the RI scheme, 
however, the specific choice of eq. floTf ) has the advantage that the vector and 
axial-vector currents, renormalized according to the same rules used for the 
four-fermion operators, satisfy automatically the relevant Ward identities. 
This is true for all the regularizations used in the intermediate steps and 
thus provides a useful check of the calculations. The validity of the Ward 
identities among renormalized quantities is not a-priori guaranteed, and it 
does not occur, for instance, in the HV- or DRED-MS schemes, because of 
the chiral symmetry breaking induced by the regularization. In the latter 
cases, the finite one- loop coefficient, entering the forward matrix element 
of the axial- vector current, does not vanish. These finite corrections are 
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compensated, in the evolution equation of the Wilson coefficients, by a term 
appearing in the two-loop current anomalous dimension. 

In order to obtain the anomalous dimension in the MS scheme, one can 
proceed in two ways. The first was explained in refs. 0-10 and the details 
will not be given here. It is based on the relation 



'MS 



B 1 - 



M) (a ) - \ (A 1 ) (lo) - PoAo - (3° X X 



(61) 



which can be derived from eqs. (155]) and fl56D by using f^jg = A . In eq. (|6l|) , 
we denoted as Aj the matrix elements restricted to the operators of the 
four-dimensional basis, and as Ai those connecting the operators of the four- 



7^ is obviously the restricted 



dimensional basis with the effervescent ones. 

matrix. The second method to obtain 7^, is by using the relation 

a J 

7$ =7^-2 (A 1 A - A A 1 ) - 2(3 A - 2/3°A^ 



(62) 



which follows from eqs. (|56|) and (|61| ) and it is also valid diagram by diagram. 
Note that the change of scheme in eq. fl62|) is equivalent to the substitution 



AlS 



■hi ~ r- 



MS 



(63) 



discussed in eq. (p9[) of subsec. |3T2. 



4.3 Checks of the calculation 

In this subsection, we illustrate several checks that have been made in order 
to verify the correctness of our calculations: 

• in the NDR-MS renormalization scheme, the ADM of the operators Qf, 
Q2 and Qf can be extracted from the results of refs. [II], 13 . They 
agree with the results presented in this paper; 

• we have computed the anomalous dimension matrix both in MS and in 
RI and verified that the result satisfy eq. (|55D , i.e. that we get exactly 
the same G in the two schemes; 
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• by computing the two-loop diagrams both in the MS and in the FRI 
schemes, we verified that the relation in eq. ( |62D holds, as we mentioned 
above, diagram by diagram [O . 



5 The anomalous dimension matrix 



In this section, we give the results for the anomalous dimension matrix in 



the Feynman-gauge RI scheme, which will be defined precisely in subsec. [5T 



and the matrices necessary to pass from FRI to a) NDR-MS, as defined 
in refs. f7[|- ||13|| ; b) the Landau-gauge RI scheme, which is the most suit- 
able for the calculation of the matrix elements on the lattice, using operators 
renormalized non-perturbatively fl5|-||17|. The FRI scheme is presented only 



because it is the simplest for doing the perturbative calculations. In practical 
cases, we expect that the MS and LRI schemes will be used for phenomeno- 
logical applications. 



5.1 The anomalous dimension matrix at LO and at the 
NLO in FRI 

In this subsection, we present the results of the leading order ADM and of 
the next-to-leading order ADM in the FRI scheme. 

The results in FRI have been obtained by computing the one- and two- 
loop Feynman diagrams shown in figs. 0-0 in the Feynman gauge. At one 
loop, we have taken the external quark momenta as indicated in fig. [IJ in 
the two-loop case, when the external subdiagrams are D\, D2 or D%, the 
external momenta have been chosen as the corresponding ones in fig. [1|. The 
field renormalization constant is computed in FRI according to eq. (|60|). 

Although the results in RI only depend on the external momenta and the 
gauge, but not on the regularization, we specify that we did the calculation 
using NDR. The choice of the external momenta may appear rather strange, 
since it is different for the different one-loop diagrams. It is, however, par- 
ticularly convenient for the perturbative calculation. In the MS scheme, the 
results for the ADM are not affected, since the independence of the external 
states is valid diagram-by-diagram. In tables [l] and § we give the complete 
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Figure 1: One-loop Feynman diagrams. We show the external quark momenta 
chosen to obtain the results in the FRI scheme. In the LRI scheme, all the 
external momenta are equal to p. 



list of the single poles necessary to compute the one- and two-loop anomalous 
dimension matrix. 

In order to present results for the ADM, we expand the coefficients of 
Ifrii written as in eq. ([14]) , in powers of a s 



4tt 1 (4tt) 
4tt (4tt) 



B = TZ^i + 773^2 + • • • (64) 



Aix 1 (Air) 2 2 

and give the expression for these quantities. 

At one-loop the ADM is independent of renormalization scheme, external 
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Table 1: Single pole contributions of the one- and two-loop diagrams to the 
ADM in the FRI scheme. In this scheme the double poles are absent. The 
label i — > j denotes the indices of the mixing-matrix. Thus 4 — > 5 corresponds 
to [Z\ ]45 ; i.e. the mixing of the bare operator with Q$, without colour 
factors. The first column refers to the diagram labels defined in fig. [7] and [| ; 
the second column to the diagram multiplicity. For the last three diagrams, 
we indicate separately the term proportional to N c or nj coming from the 
gluon vacuum-polarization of the internal gluon line, see fig. [|. 
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Table 2: Single pole contributions to the ADM of the one- and two-loop self- 
energy diagrams in the FRI scheme. For S 2 the contributions proportional 
to N c (3rd column) and rif (4th column) are shown separately. 
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Figure 3: Two-loop Feynman diagrams. 
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states and gauge. The results in this case are the following 
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Our results agree with those obtained in ref. 0. 

For the two-loop ADM we obtained (with /3° computed from eq. ( f4"0| ) 
A = l) 
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For Jfju , we have shown explicitely those terms, proportional to which 
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5.2 Relation between FRI and other renormalization 
schemes 

In this subsection, we give the recipe to pass from FRI to other schemes 
which may be useful for practical applications: LRI, for lattice calculations, 
and the standard MS NDR scheme. All we need to know is the shift matrix 

JlRI = JfRI + ^FRI 

Jms = JFRi+f T FRI -r?—. (69) 

From the knowledge of J in a given renormalization scheme, we can imme- 
diately obtain the evolution matrix W[pi, M] using eqs. (^) and (|23|). 

As discussed in subsec. [3~2] , the renormalization scheme is completely 
defined by the matrix f of eq. (^), computed for given external momenta 
and gauge. We choose quarks with equal momentum p as external states and 
the Landau gauge. The field renormalization constant is computed in LRI 
according to eq. (|60"D , which gives Z^ — l. We denote as r FRI , f LRI and f^g 
the three cases considered here: f l R i is obviously zero and fj^g is a 20 x 20 
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matrix because this regularization (and consequently the corresponding MS- 
renormalization scheme) does not respect chiral and Fierz symmetries. By 
denoting with f ++ the 5x5 sub-matrix in the Qf sector and similarly for 
f + ~, r~ + and f , we get (see also ref. [|14j] for the operators Q^ 2)3 ), 
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